The SMAA-PROMETHEE methods 
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Abstract: PROMETHEE methods are widely used in Multiple Criteria Decision Aiding (MCDA) 
to deal with real decision making problems. A crucial aspect of the classical PROMETHEE 
methods is the choice of criteria weights. In this paper, we propose to apply the Stochastic 
Multiobjective Acceptability Analysis (SMAA) to the classical PROMETHEE methods and to 
the bipolar PROMETHEE methods in order to explore the whole set of weights compatible with 
some preference information provided by the Decision Maker (DM). A didactic example describes 
the application of the presented methodology to a student evaluation problem. 
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1 Introduction 



PROMETHEE is a well-known family of Multiple Criteria Decision Aiding (MCDA) outranking methods 
(see [13] for a survey on MCDA and [H [5] for a survey on PROMETHEE methods) and widely used in order 
to deal with several types of real world problems [3]. The basic assumption of the PROMETHEE methods, 

^ | the most known of which are PROMETHEE I and II, is that the set of criteria is mutually preferentially 

H ■ 

C?.' independent [42J. In several decision making problems, this assumption is in general violated since criteria 
are frequently interacting with each other, which means that it is possible to observe a certain degree of 
redundancy or synergy among them. For example, when evaluating sport cars, the criteria maximum speed 
and acceleration may be considered redundant because, in general, speed cars have also a good acceleration. 
Therefore, even if these two criteria can be very important for a DM liking sport cars, their comprehensive 
importance is smaller than the sum of importance of the two criteria considered separately. From the 
other side, criteria maximum speed and price lead to a synergy effect because a speed car having also a 
low price is very well appreciated. For such a reason, the comprehensive importance of these two criteria 
should be greater than the sum of importance of the two criteria considered separately. In these two cases 
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the aggregation of the preferences is done by using non-additive integrals the most know of which are the 
Choquet integral [Tj and the Sugeno integral [37] (see [TTJ [2TJ 12Q] for the application of non-additive integrals 
in MCDA). 

In other cases, the importance of a criterion expressing preference in favour of an alternative can be reduced 
by the presence of another criterion opposing to it. For instance, when considering again the car evaluation 
problem, the weight of criterion speed in favor of an alternative a can be reduced by the presence of 
criterion aesthetic opposing to it. In these cases, the generalization of Choquet and Sugeno integrals, that 
are the bipolar Choquet integral and the bipolar Sugeno integral can be used to aggregate the preferences 

pa QS EH Egg. 

In [9] and [8], the authors presented the bipolar PROMETHEE methods extending the PROMETHEE 
I and II methods to the case in which criteria are not independent, proposing to use the bipolar Choquet 
integral to aggregate the preferences. 

In order to apply these methods, the Decision Maker (DM) has to provide inter-criteria and intra- 
criterion preference information on the parameters involved. This type of preference information can be 
given directly or indirectly. Direct preference information consists of asking the DM to provide precise 
values of the parameters, while indirect preference information consists of asking the DM to provide some 
comparisons regarding reference alternatives from which the values of the parameters should be elicited. 
The indirect preference information demands less cognitive effort from the DM and, for this reason, it is 
widely used in outranking methods (see for example [TO], [3H [351 [231 [29]). 

Recommendations provided by the PROMETHEE methods are dependent on the values given to the 
considered parameters; in fact, generally, different sets of parameters lead to different comparisons between 
the alternatives. For this reason, many methods have been proposed in literature to elicit sets of parameters 
compatible with some preference information provided by the DM or to get values of the parameters giving 
to an alternative the best rank. In the following we cite only some of these contributions. 
Considering an MCDA additive method, and getting consequently a ranking of the considered alternatives, 
Mareschal |32] makes a sensitivity analysis to obtain stability intervals for the criteria weights such that the 
method provides the same ranking whichever the choice of the weights inside these intervals. Sun and Han 
|38j solve a linear programming problem to find the most discriminant set of weights compatible with the 
preference information provided by the DM. Wolters and Mareschal |43] propose a method to study how the 
ranking of the alternatives is sensitive to the changes on the weights and on the alternatives' evaluations, 
and which modifications on the importance weights are necessary to make an alternative the best. Eppe and 
De Smet [11] study the weights that best represent the preferences of the DM and the type of information 
structure giving the ranking closest to the ranking provided by the DM Eppe et al |12j . use a meta-heuristic 



in order to induce parameters for the PROMETHEE II model starting from some preference information 
provided by the DM. 

If the DM decides to give indirect preferences, Robust Ordinal Regression (ROR) [251 El E3 [23J can be 
applied to explore the whole space of parameters compatible with the preference information provided by 
the DM. ROR is a family of MCDA methods taking into account simultaneously all the sets of parameters 
compatible with some preference information provided by the DM considering two preference relations, one 
possible and one necessary. The necessary preference relation holds between two alternatives a, b if a is at 
least as good as b for all sets of parameters compatible with the information provided by the DM, while 
the possible preference relation between two alternatives holds if a is at least as good as b for at least one 
set of parameters compatible with the preference information provided by the DM. ROR for preferences 
with interacting criteria modelled through the Choquet integral has been considered in [2j. ROR has been 
applied also to the classical PROMETHEE methods in [23] and to the bipolar PROMETHEE methods in 
[8j. 

Similarly to the ROR, the Stochastic Multiobjective Acceptability Analysis (SMAA) [301 EI] explores 
the whole set of parameters compatible with some preference information provided by the DM. SMAA is 
a family of MCDA methods taking into account imprecision or lack of the data considering probability 
distributions on the space of criteria weights and on the space of alternatives' evaluations. SMAA has been 
applied to preferences with interacting criteria modelled through the Choquet integral in pQ. 
In this paper we propose an integrated approach between SMAA and the classical and the bipolar PROMETHEE 
methods giving arise to the SMAA-PROMETHEE methods. From one side, integrating SMAA and PROMETHEE 
I methods permits to compute the frequency of the preference, indifference or incomparability between two 
alternatives From the other side, integrating SMAA and PROMETHEE II methods permits to study the 
different positions an alternative can get when varying the parameters of the model. 

The paper is organized as follows. In the next Section we shall describe the classical PROMETHEE 
I and II methods while in Section [3J we recall the basic concepts of the bipolar PROMETHEE methods; 
Section [5] contains the description of SMAA methods when the preference model is an additive function; in 
Section [5] we shall present the SMAA-PROMETHEE methods and we shall describe its link with the ROR: 
Section [6] contains a didactic example while some conclusions and future directions of research are presented 
in the last section. 



2 The classical PROMETHEE methods 



Let us consider a set of actions or alternatives A evaluated with respect to a set of criteria G = {gi, . . . ,g n }, 
where gj : A ->• R, j G J = {1, . . . , n} and |A| = m. PROMETHEE gj E] is a well-known family of MCDA 
methods that aggregate preference information of a DM through an outranking relation. Considering for 
each criterion g~ a weight Wj (representing the importance of criterion gj within the family of criteria G) , an 
indifference threshold qj (being the largest difference dj(a, b) = gj{a) — gj(b) compatible with the indifference 
between alternatives a and b on criterion gj), and a preference threshold pj (being the minimum difference 
dj(a,b) compatible with the preference of a over b on criterion gj), PROMETHEE methods build a non- 
decreasing function of dj(a,b), whose formulation (see [I] for alternative formulations) can be stated as 
follows 



Pj(a,b) = < 



if dj(a, b) < qj, 



if Qj < dj(a,b) < Pj , (1) 



dj (a,b)—qj 
Pj-1j 

1 if dj(a, b) > pj. 

Pj(a, b) represents the degree of preference of a over b on criterion gj. 

For each ordered pair of alternatives (a, b) € A x A, PROMETHEE methods compute 

= ^2 w 3 p j( a > b ) 

representing how much alternative a is preferred to alternative b taking into account the whole set of criteria. 
ir(a,b) can assume values between and 1 and obviously the greater the value of 7r(a, 6), the greater the 
preference of a over b. 

In order to compare an alternative a with all the other alternatives of the set A, PROMETHEE methods 
consider the negative and the positive flow of a defined in the following way: 

<fr(a) = — — - y~] n(b,a) and 0+(a) = — l — S~] vr(a,6). 
m — 1 ^-^ m — 1 ^-^ 

b£A\{a} beA\{a} 

<f>~(a) represents how much the alternatives of A \ {a} are preferred to a; the lower (ft~(a) the better al- 
ternative a is; (j) + (a) represents how much a is preferred to the alternatives of A \ {a}; the greater (j) + (a), 
the better a is. PROMETHEE II computes also the net flow <j>(a) = <j) + {a) — <fi~(a) for each alternative 
a £ A. Taking into account the negative and the positive flows, PROMETHEE I builds a preference (V 1 ), 
an indifference (I 1 ) and an incomparability (R- 1 ) relation on the set of alternatives A, while considering 
the net flow, PROMETHEE II builds a preference {V 11 ) and an indifference (X n ) relation on the set of 



alternatives A: 

aV 1 b iff <£ + (a) > <& + (&), &~{a) < &~{b) and at least one of the two inequalities is strict, 
< aX 7 6 iff $+(a) = $+(6) and $-(a) = 
aRJb otherwise. 

aP^b iff $(a) > $(6), 

< 

aX /7 6 iff $(a) = $(6). 

PROMETHEE I and PROMETHEE II provide, respectively, a partial and a complete order of the alterna- 
tives in A. 

3 The bipolar PROMETHEE methods 

In order to take into account the bipolar preferences, the bipolar PROMETHEE methods [HJ E] build for 
each criterion j G JT, and for each pair of alternatives (a, 6) S i x i, the bipolar preference function 

f Pj(a,b) if Pi(a,&) > 
Pf(a,6) = A,> (2) 

[ -Pj(b,a) if Pj(a,6) = 0. 

It is representative of the preference of a over b on criterion gj and it is straightforward proving that 
Pj(a, b) = —Pj(b, a) for all j G J7" and for all a,i»£ A 

In the bipolar PROMETHEE methods, the aggregation of the preference functions Pf(a, b) is done by using 
the bipolar Choquet integral. 

The bipolar Choquet integral is based on a bicapacity [HI [19], being a function ft defined on P(J) = 
{(C,D) : C,D C J and C n D = 0} and having value on [-1, 1] such that: 

• pL{$,J) = -l,p,(J,<ft) = 1, £(0,0) = (boundary constraints), 

• £(C, -D) < fi(E,F) for all (C,D),(E,F) G P(<7), such that C C E and D D F (monotonicity 
constraints) . 

According to [22J, in the bipolar PROMETHEE methods we consider a decomposed bicapacity £ [23] such 
that £(C,P) = /x+(C,P) - iT{C,D) for all (C,D) G P{J) where : P(J) -»■ [-1,1] satisfy the 

following conditions: 

• ^(J-, 0) = 1 and /x+(0, P) = 0, VP C J, 



• M -(0, J) = 1 and a*" (-8,0) =0, VBC J, 

• H + {C,D) < H + {E,F) for all F) G :CCE, D D F, 

• fi~(C,D) < fi~(E,F) for all (C,D), (E,F) £ P{J) :C^E, DQF. 

The interpretation of the functions fi + and jjT is the following. Given the pair (a, b) £ A x A, let us 
consider (C, -D) € P{J) where C is the set of criteria expressing a preference of a over 6 and D the set of 
criteria expressing a preference of 6 over a. In this situation, fi + (C,D) represents the importance of criteria 
from C when criteria from D are opposing them, and fi~ (C, D) represents the importance of criteria from 
D opposing C. Consequently, fi(C, D) represents the balance of the importance of C supporting a and D 
supporting b. 

After reordering in a non-decreasing way the preference functions for each criterion j € J ', (\P%Ja, b)\ < 
\P^(a, 6)| < ... < \P B ^(a, 6)| < ... < \P^(a, b)\), on the basis of the decomposed bicapacity fi = — \i~ , 
we can give the following definitions of the bipolar comprehensive, positive and negative preference: 



Ch B (P B (a,b),jl)= \P§)(aM[KC {3) ,D {3) ) - ttC {3+1)l D {]+1) ) 



(3) 



Ch B+ (P B (a,b),fi + ) = \P { %(aM[^(C U) ,D U) ) - ^(C ij+1) ,D U+1) ) 



(4) 



Ch B -(P B (a,b),tJL-) = \P§){aM[^(C (3) ,D (3) ) - ^{C {]+l) ,D (3+l) ) 

jej> 



(5) 



where P B (a,b) = [P B (a,b), j € j\, J> = {j € J : |P ( f } (a,6)| > 0}, C {j) = {i € J> : P B (a, 6) > 
|P ( f } (a,6)|}, D {J) = {zej> : -Pf(a,b) > |P^(a,6)|} and C (n+1) = D (n+1) = 0. 



Ch B (P B (a,b), ft) gives the comprehensive preference of a over 6 and it is equivalent to ir(a, b) — ir(b,a) 
in the classical PROMETHEE method. It is, therefore, reasonable expecting that Ch B (P B (a, 6), (t) = 
-Ch B (P B (b, a), p) for all a, b € A. 

Ch B+ (P B (a,b), ) and Ch B ~ (P B (a,b), give, instead, how much a outranks b (considering the reasons 
in favor of a) and how much a is outranked by b (considering the reasons against a). Also in this case, it is 
reasonable expecting that Ch B+ (P B (a,b), = Ch B ~ (P B (b, a), for all a, b € A. 
From the definitions above, it is straightforward proving that, for all a, b <G A, 



Ch B (P B (a,b),fL) = Ch B+ (P B (a,b), f i + )-Ch B -(P B ( y a 1 b),n-) 



(6) 



Similarly to the classical PROMETHEE methods, also in the bipolar PROMETHEE methods one can 
define the bipolar positive, negative and net flows for all a € A: 



B+ («) = 7 £ Ch B+ (P B (a,b),^) (7) 

bGA\{a} 

cp B -(a) = ^— Y, Ch B -(P B (a,b),^) (8) 
m — 1 

bGA\{a} 

<t> B {a) = -^— Y Ch B (P B (a,b),fL) (9) 
m — 1 

b£A\{a} 

and by equation ([6]), it follows that 4> B {a) = <fi B+ {a) — (j) B ~(a) for all a £ A. 

Equivalently to PROMETHEE I method, the bipolar PROMETHEE I method builds a preference (V^), an 
indifference and an incomparability (7£|j) relation, while equivalently to the classical PROMETHEE 

II method, the bipolar PROMETHEE II method builds a preference {VU) and an indifference (Ig) relation 
on the set of alternatives A: 

aVgb iff <S> B+ (a) > 3> s+ (6), <S> B -(a) < <5> B -(b) and at least one of the two inequalities is strict, 
< algb iff <5> B+ (a) = $> B+ {b) and $ B -(a) = 
aR/g b otherwise. 

aVgb iff <S> B (a) > 

< 

aXgb iff <5> B {a) = <5> B (b). 

Because the use of the bipolar Choquet integral, and therefore of the bipolar PROMETHEE methods, 
is based on a bicapacity which assigns numerical values to each element of P(J), in [8] we have considered 
only the 2-additive bicapacities [18^ I16j . being a particular class of bicapacities. We give the following 
decomposition of the functions [i + and [i~ previously defined: 

. / z + (C, J D) = ^a + ({i},0)+ Y, + Y ° + (0'M*}) 

j&C {j,k}QC jeC, k&D 

. /,-(C, J D) = 5>-(0,{j})+ Y *~(W<}) + Y »-({#»{*}) 

jeD {j,k}CD jeC, k&D 

where: 

• a + ({j},0) represents the power of criterion gj by itself expressing positive preferences; this value is 
always non negative; 



• a + ({j, k}, 0) represents the interaction between gj and gk when they are in favor of the preference of 
a over b; this value is zero if there is no interaction; it is positive (negative) when there is a synergy 
(redundancy) effect when putting together gj and 

• a + ({j}, {k}), represents the power of criterion against criterion gj, when criterion gj is in favor of 
a over b and g^ is against to the preference of a over b; this leads always to a reduction or no effect on 
the value of since this value is always non-positive. 

An analogous interpretation can be applied to the values a~(0, {j}), a~(0, {j, k}), and a~({j}, {k}). 

Let us remember that positive and negative interactions, and power of one criterion against another 
correspond to the mutual strengthening effect, the mutual weakening effect and the antagonostic effect 
considered for the ELECTRE methods in [H] . 

From [i + and fj,~ we get the following 2-additive decomposable bicapacity fx 

A(tf,s) = S>?-£«7+ E 4- E a J k + E 4*- E V (io) 

jec jeD {i.fcjcc {j,A;}cD jec*, fceD jec*, feeD 

where at, a+ fc , a+ fc denote a+({j},0), a + ({j*,fc},0), a+({j}, {A;}), and a~, aj k , aTj fe denote a~(0,{j}), 
a~(0,{j,fc}), a"({j},{/c}), respectively. 

On the basis of a 2-additive decomposable bicapacity, the monotonicity and the boundary conditions 
previously considered can be rewritten as follows (see [8]): 

• 4 + E 4 + E > °> v 3 e V(C u {j}, D) e P(^), 

fcGC keD 

• E °Jj ^ °> V 3 G ^ V(C, D U {j}) G P(J), 

• °7 + E a 7* + E a %- ^ °> v J G ^ v ( c ' D u i3}) e 

• E a ^^ ' VjGj, V(CU{j},D)GP(J), 

• E4+ E 4^ = 1' 

ie^ {j,k}cj 

• E a 7+ E 7fc = 1 - 

Considering a 2-additive decomposable bicapacity, in [8] we have given the following theorems defining the 
comprehensive bipolar positive flow and the comprehensive bipolar negative flow; before we have provided 



conditions so that Ch B \P B (a,b), p.) = -Ch B (P B (b,a), fi) and Ch B+ (P B (a, b), fx+) = Ch B - (P B (a,b), 
for all a, b G A. 

Theorem 3.1. (see [8]) Given a 2-additive decomposable bicapacity ft, then for all x G W 1 

1. Ch B+ (x,fi + ) = ^ af k mxn{xj,x k } + a^ k min{xj, -x k } 

j£j,Xj>0 j,keJ,j^k,Xj,x k >0 j,k£j,j^k,Xj>0,x k <0 

2. Ch B ~(x,fi~) = - a j x j~ ^ aj k max{xj,x k } - ^ a~j fc max{-Xj, x k } 

j£j,Xj<0 j,keJ,j^k,Xj,x k <0 j,k£j,j^k,Xj>0,x k <0 

Theorem 3.2. (see Given a 2-additive decomposable bicapacity fi, Ch B (P B (a,b), p,) = —Ch B (P B (b,a),fi) 
and Ch B+ (P B {a, b),fi + ) = Ch B -(P B (b, a),fT) for all a,b G A iff 

1. for each j G J, = a~ , 

2. for each {j, k} C J, a+ fc = aj k , 

3. for each j, k G J, j ^ k, a1j fe = a~^. 

4 SMAA 

Chosen the decision model and assigned the values of its parameters, it is straightforward obtaining some 
recommendation regarding the problem at hand. In general, the DM is not able to provide directly these 
parameters or she is not worth of their meaning; besides, in real world decision making problems, data 
regarding the alternatives' evaluations can be missing or imprecise. 

SMAA is a family of MCDA methods taking into account imprecision or lack of data in the problem at 
hand (see |39] for a comprehensive survey on the use of SMAA methods in MCDA) considering a probability 
distribution fw over the space of all compatible weights W and a probability distribution f x over the space 
X C ]J mxn of the alternatives' evaluations. 

Without any preference information provided by the DM, the weight space is defined as follows: 

< w G R n : Wj > 0, Vj G J, and = l\ . 

jej J 

If the DM is able to provide some preference information, then the space W is restricted by the constraints 
translating these preference information. Consequently, in the following, with W we shall denote the set of 
weights compatible with the preference information provided by the DM. Now let us describe some SMAA 
concepts considering the additive value function, u(xi,w) = '^2 w j9j( x i) with X{ G A and w G W, as 
preference model. 



Because u(x{,w) provides a complete ranking of the alternatives, for each £ in \ an d w in W, SMAA 
computes the position reached by alternative Xi € A as 

rank(i,£,w) = 1 + ^p(u(£ fc ,w) > u(&,w)), 

where p(true) = 1 and p(false) = 0. Besides, for each £ 6 x, SMAA computes the favourable rank weights 
of alternative x% € A 

Wf(£) = 6 ^ : rank(i,£, w) = r} 

being the set of possible weights giving to alternative X{ the position r = 1, . . . , m in the final rank. 

On the basis of the favourable rank weights, SMAA computes the rank acceptability index and the central 

weight vector. 

• The rank acceptability index 

b r i = [ fx(0 I fwH dw d£ 

describes the share of parameters giving to alternative Xi the position r in the obtained final ranking; in 
particular, b\ measures the variety of parameters making alternative Xi the most preferred. Obviously, 
the best alternatives are those having rank acceptability index greater than zero for the first positions 
and rank acceptability index close to zero for the lower positions. 

• The central weight vector 

w i= tt [ fx(0 I fw(w)w dw di 

describes the preferences of a typical DM that makes alternative Xj the most preferred. 

5 SMAA-PROMETHEE 

Even if SMAA has already been applied to outranking methods and in particular to ELECTRE methods 
[281 HQ], to the best of our knowledge, no attempt of applying the SMAA methods to PROMETHEE has 
been done in literature before. In [8], we pointed out that the bipolar PROMETHEE methods are reduced 
to the classical PROMETHEE methods if = a^j fc = aTj fc = for all j, k £ J~,j ^ k and in this case 
the parameters dj in the bipolar PROMETHEE methods are the weights Wj of the classical PROMETHEE 



methods. For this reason, without loss of generality, in this section we shall apply SMAA only to the bipolar 
PROMETHEE methods. 

We shall suppose that parameters aj, a1j fe , and a~^ k , are unknown for all j,k G J and, at the same time, 
we shall assume that indifference and the preference thresholds as well as the alternatives' evaluations are 
already known. 

In the bipolar PROMETHEE methods, we consider the following information provided by the DM and 
their representation in terms of linear constraints: 

1. Cora-paring pairs of actions locally or globally. The constraints represent some pairwise comparisons on 
a set of training actions. Given two actions a and b, the DM may prefer a to b, b to a or be indifferent 
to both: 

(a) the linear constraint associated with aVb (a is locally preferred to b) is: 

Ch B (P B (a,b),fi) > 0; 

(b) the linear constraints associated with aV l B b (a is preferred to b with respect to the bipolar 
PROMETHEE I method) are: 

$ B +(a) > $ B +(6), 
§ B -(a) < $ B ~{b), 

$ B+ (a) - $ B -(a) > $ B+ (b) - $ B -(b), 

(c) the linear constraint associated with aV^b (a is preferred to b with respect to the bipolar 
PROMETHEE II method) is: 

$ B (a) > $ B {b) 

(d) the linear constraint associated with alb (a is locally indifferent to b) is: 

Ch B (P B (a,b),fi) = 

(e) the linear constraints associated with aI l B b (a is indifferent to b with respect to the bipolar 
PROMETHEE I method) are: 

$^+(a) =$ B +(b), 1 
$ B -(a) = $ B -(6), J 



(f) the linear constraint associated with aXgb (a is indifferent to b with respect to the bipolar 
PROMETHEE II method) is: 

$ B ( a ) = q> B (b) 

2. Comparison of the intensity of preferences between pairs of actions. The constraints represent some 
pairwise comparisons between pairs of alternatives on a set of training actions. Given four actions a, 
b, c and d: 

(a) the linear constraints associated with (a, b)V(c, d) (the local preference of a over b is larger than 
the local preference of c over d) is: 

Ch B (P B (a,b),fi) > Ch B (P B (c,d),fi) 

(b) the linear constraints associated with (a, b)T(c, d) (the local preference of a over b is the same of 
local preference of c over d) is: 

Ch B (P B (a,b),fi) = Ch B (P B (c, d),jl) 

Importance of criteria. A partial ranking over the set of criteria J may be provided by the DM: 

(a) criterion gj is more important than criterion g^, which leads to the constraint aj > a^; 

(b) criterion gj is equally important to criterion g^, which leads to the constraint aj = a&. 

The sign of interactions. The DM may be able, for certain cases, to provide the sign of some inter- 
actions. For example, if there is a synergy effect when criterion gj interacts with criterion 5%, the 
following constraint should be added to the model: ajk > 0. 

Interaction between pairs of criteria. The DM can provide some information about interaction between 
criteria: 

a) if the DM feels that interaction between gj and gk is greater than the interaction between g p and 
g q , the constraint should be defined as follows: \ajk\ > \a pq \ where in particular: 

• if both couples of criteria are synergic then: ajk > a pq , 

• if both couples of criteria are redundant then: ajk < a pq , 

• if (j, k) is a couple of synergic criteria and (p, q) is a couple of redundant criteria, then: 



• if (j, k) is a couple of redundant criteria and (p, q) is a couple of synergic criteria, then: 

b) if the DM feels that the strength of the interaction between gj and is the same of the strength 
of the interaction between g p and g q , the constraint will be the following: \ajk\ = \a pq \ and in 
particular: 

• if both couples of criteria are synergic or redundant then: ajk = a pq , 

• if one couple of criteria is synergic and the other is redundant then: ajk = — a pq , 

6. The power of the opposing criteria. Concerning the power of the opposing criteria several situations 
may occur. For example: 

a) when the opposing power of <j% is larger than the opposing power of gh, with respect to gj, 
which expresses a positive preference, we can define the following constraint: a1j fc < a1j h (because 
at h < and aj^ h < for all j, k with j / k); 

b) if the opposing power of expressing negative preferences, is larger with gj rather than with 
gh, the constraint will be a^ k < a^ k . 

The constraint translating the preference information of the DM along with the symmetry, boundary 
and monotonicity conditions make the following set of constraints in which strict inequalities have been 
transformed in weak inequalities by adding an auxiliary variable e: 



if aV^b 



Ch B (P B (a,b),jl) > e if aVb, 
1> B +(a) > $> B +(b), 
* s -(a) < <S> B ~(b), 

$ B +(a) - <!> s -(a) > $ s +(6) - $ s -(6) + s 
<f> B (a) > <S> B (b) + e if aV'Jb 

Ch B {P B (a,b),fl)>Ch B {P B (c,d),fi) + e if (a,b)V(c,d), 
cij — a,k > e if 

l a jfcl — \ a pq\ > £ if {i, g}, (see point 5. a) above) 

l a jfcl = | otpg | if {j,k}X{p,q}, (see point 5.b) above) 
a jk > £ if there is synergy between criteria j and k, 
®jk < — £ if there is redundancy between criteria j and k, 
djk = if criteria j and k are not interacting, 
Power of the opposing criteria of the type [6] 
a~h, — a + , > e, 
aft, — a + , , > e, 

]\k p\k — 

Symmetry conditions (Theorem l3.2l l: 



Ch B (P B (a,b),fi) = if alb, 



$ s +(a) = <S> B +(b), 
$B-(a) = <3> B ~(b) 



if algb 



<J> s (a) = $ s (b) if aZj/fe 

Ch B (P B (a,b),fi) = Ch B (P B (c,d),ft) if (a, b)X(c,d), 
a>j = ife if i^fc, 



> E, 



-a |, > e, 

p\k — 



4* 



:a" u , Vj,kej,j^k 



Boundary and monotonicity conditions: 

jej {j,k}^J 
a-j > V j e J, 

a 3 + J2^k+J2 a l\k ^ °> VjeJ, V(C U {j}, D) G P(J), 
kec feso 

aik + E a hH ^ °. v J e ^ D u 0'}) e P( J). 



' / * Jk 
keD 



In order to progress from the simplest to the most sophisticated method, one could check firstly if the 
preferences of the DM are compatible with the classical PROMETHEE methods; if this is not true, then 
one can pass to the bipolar PROMETHEE methods. 

Remark 5.1. The DM could start directly by checking if the bipolar PROMETHEE methods can restore 
her preference information without looking at the compatibility of the classical PROMETHEE methods; in 
the same way, she can check if the bipolar PROMETHEE methods can restore her preferences even if the 
classical PROMETHEE methods are already able to restore them. 

In order to check if the preference information provided by the DM can be restored by the classical 
PROMETHEE methods, one has to solve the following optimization problem: 



max e = e\ 
E B 



(ii) 



If E B is feasible and e\ > 0, then the classical PROMETHEE methods are able to restore the preference 
information provided by the DM, otherwise, solve the following optimization problem 

max e = e\ s.t. 

(12) 

E B 

If E B is feasible and e\ > 0, then the bipolar PROMETHEE methods can restore the preference information 
provided by the DM, otherwise one can try to discover which constraints cause the incompatibility by using 
some technique presented in [33J. 

Since the constraints in E B are linear, we propose to apply the Hit- And- Run method |36|, |4"T] to sample 
sets of parameters compatible with the preference information provided by the DM. 

The Hit-And-Run sampling begins with the choice of one point inside the polytope delimited by the con- 
straints translating the preference information provided by the DM. At each iteration, it is sampled from 
the unit hypersphere a random direction that, with the considered position, generate a line. Finally, the 
sampling of one point inside the segment whose extremes are the intersection of the line with the bounds 
and the current point is done. 

Proceeding in this way, we sample a given number of sets of parameters and, applying the PROMETHEE 
methods considering each time one of these compatible sets of parameters, we store: 

1. preference Vg, indifference Zg and incomparability / Rl B relations obtained by applying the bipolar 
PROMETHEE I method, 

2. preference Vg and indifference Xjj relations obtained by applying the bipolar PROMETHEE II 
method, 

3. ranking of the alternatives by applying the bipolar PROMETHEE II method, 

4. set of parameters sampled in the iteration. 

Sampled the fixed number of sets of parameters, we get the following results: 

• Considering the ranking obtained by the bipolar PROMETHEE II method for each parameters' sam- 
pling, one can compute the rank acceptability index K, the preference matrix P and the indifference 



matrix /. In particular, b\ describes the share of parameter values giving to alternative x\ rank r; 
the element pij of the preference matrix P is the frequency of the preference of alternative Xi over Xj 
while the element of the indifference matrix / in the position is the frequency of the indifference 
between Xi and Xj. Besides, one can compute the central weight vector w\ and the barycenter w c 
representing respectively the typical preferences allowing to alternative Xi to get the first rank and the 
average preferences of the DM; 

Remark 5.2. We observe that wf and w c are vectors of R™ if we apply SMAA to the classical 

i 2 

on —n 

PROMETHEE methods, while they belong to R~ if we apply SMAA to the bipolar PROMETHEE 
methods ^ n2 ~ n _ n _|_ @ 2 _|_ D n % where n are the weights of criteria, C n> 2 is the number of the 
parameters regarding the interaction between couples of criteria and P> n ,i * s the number of parameters 
of the type a ^ k - We did not consider the parameters aj, k because a1j fc = o,^.). 

• Considering the bipolar PROMETHEE I method, for each couple of alternatives Xi,xj £ A one can 
compute the frequency of the preference of x% over Xj , the frequency of the preference of Xj over x% or 
still the frequency of the indifference or incomparability between Xi and Xj) 

In [SJ, the authors have applied the ROR [27] to the bipolar PROMETHEE methods. ROR is a family 
of MCDM methods tacking into account simultaneously all the sets of parameters compatible with prefer- 
ence information provided by the DM building a necessary and a possible preference relation. Given two 
alternatives Xi and Xj, we say that xi is necessarily preferred to Xj, and we write xi ^ N Xj, if xi is at least 
as good as Xj for all compatible sets of parameters, while Xi is possibly preferred to Xj, and we have write 
X{ ^ p xj, if Xi is at least as good as Xj for at least one compatible set of parameters. ROR has been already 
applied to both classical PROMETHEE methods [23] and to the bipolar PROMETHEE methods [8]. 
Because the computation of the necessary and possible preference relations involve to solve two linear pro- 
gramming problems for each pair of alternatives in A, we could use the SMAA-PROMETHEE methods to 
approximate the two preference relations. In fact: 

• if Xi is necessarily preferred to Xj, then the sum of the frequency of the preference of Xi over Xj and 
of the frequency of the indifference between Xi and Xj is 100%, 

• if the frequency of the preference of Xi over Xj is greater than zero, then Xi is possibly preferred to Xj. 

The vice versa of these statement are not true because, from one side, also if the frequency of the preference 
of Xi over Xj is the 100%, there could exists one non sampled set of compatible parameters for which Xj is 
preferred to Xf, from the other side, also if Xi is possibly preferred to xj, it is possible that for all sampled sets 



of parameters Xj is at least as good as Xi and therefore the frequency of the preference of xi over Xj is 0%. 
Observe that the larger is the sample of the sets of parameters, the better the approximation of the necessary 
and the possible preference relations obtained through SMAA, such that, in case of an enough large sample 
of parameters, one can reasonably accept the approximation of SMAA as results of ROR analysis. 

6 Illustrative example 

In this section we shall apply the SMAA-PROMETHEE method described in section [5] to the same example 
presented in [8|. 

Suppose that the Dean of a high school has to decide which student deserves to get a scholarship. Eight 
students arrived to the final selection and in order to make the "best" decision, the Dean decides to evaluate 
them with respect to three different subjects: Mathematics (M), Physics (P) and Literature (L). Their 
evaluations are given on a [0,20] scale as shown in Table [TJ The Dean has some preferences regarding 
two couples of students: S7 and S2 from one side and S5 and sq from the other side. She states that, 
she prefers sj to S2 and S5 to sq. These local preference information are translated by the constraints 
Ch B (P B (s 7 ,s 2 ),fi) > and Ch B (P B (s 5 , s 6 ), (1) > respectively. 

Table 1: Evaluations of the students 



Students 


Mathematics 


Physics 


Literature 


si 


16 


16 


16 


S2 


15 


13 


18 


■S3 


19 


18 


14 


S4 


18 


16 


15 


SB 


15 


16 


17 


S6 


13 


13 


19 


S7 


17 


19 


15 


S8 


15 


17 


16 



As explained in the previous section, at first we check if the preferences of the Dean are compatible with 
the classical PROMETHEE methods. In this case we get that the set of constraints E B ' is feasible and 
e\ > 0. Therefore, the classical PROMETHEE methods can restore the preference information provided by 
the Dean. 

Since the Dean is interested to get a general overview about the eight students, we present to her the results 
in Tables El and [3] got by applying SMAA to the classical PROMETHEE I and II. 

Sampling 100000 sets of parameters compatible with the preference information provided by the Dean, 
the results show that S3 and sj are surely the best two students. In fact, looking at the positive and the 
negative flows separately, that is to how much an alternative outranks all other alternatives and how much 



Table 2: Results obtained by applying SMAA to the classical PROMETHEE I method 



(a) Preferences in percentage (b) Incomparability in percentage 

Sl S 2 S 3 S 4 S 5 S 6 S 7 S S Si S 2 S3 S 4 S 5 S T Sg 

si 67,91 18,44 58,70 61,07 42,96 

s 2 2,09 7,46 38,99 11,85 

s 3 71.25 91,21 68,01 73,60 88,99 35,70 71,39 

s 4 63,43 81,31 65,77 73,74 27,83 55,67 

s 5 31,06 62,89 2,59 26,62 58,09 37,10 

s 6 9,52 1,75 

s 7 90.90 99,19 39,23 71,08 87,60 84,83 89,21 

s 8 47,4 62,89 1,55 38,48 62,9 60,40 



an alternative is outranked from all other alternatives, we get that students S3 is preferred to all other 
students, except to 37, with a frequency of at least 68.01%, while S7 is preferred to all other students, except 
to S3, with a frequency of at least 71.08%. Comparing S3 and S7, we get instead that S3 is preferred to S7 
with a frequency of 35.7%, S7 is preferred to S3 with a frequency of 39.23% and they result incomparable 
with a frequency of 25.07%. 

Table 3: Results obtained by applying SMAA to the classical PROMETHEE II method 



(a) Preferences in percentage (b) First rank acceptability 





si 


S2 


S3 


s 4 


S5 


se 


S7 


ss 




b i 


b ? 


b ? 


b f 


b ? 


b ? 


b i' 


b r 


si 





90,62 


4,67 


25,59 


63,43 


91,73 





47,22 


si 











38,567 


47,836 


11,913 


1,679 


0,005 


S2 


9,38 





3,40 


8,57 


15,14 


74,27 





20,61 


s 2 











5,539 


13,255 


5,815 


57,845 


17,546 


S3 


95,33 


96,6 





100 


88,78 


93,68 


53,61 


87,95 


S3 


53,612 


33,045 


3,211 


4,033 


1,164 


1,786 


3,149 





s 4 


74,41 


91,43 








70,29 


89,18 


28,62 


59,07 


s 4 





28,625 


28,896 


13,881 


2,778 


15,222 


2,468 


8,13 


S5 


36,57 


84,86 


11,22 


29,71 





100 


1,08 


37,10 


S5 


1,085 


7,143 


5,982 


18,146 


12,294 


55,35 








sg 


8,27 


25,73 


6,32 


10,82 








0,09 


8,62 


S6 





0,094 


2,794 


3,978 


1,113 


3,127 


19,802 


69,092 


s 7 


100 


100 


46,39 


71,38 


98,92 


99,91 





100 


S7 


45,303 


25,978 


28,719 

















S8 


52,78 


79,39 


12,05 


40,93 


62,9 


91,38 








ss 





5,115 


30,398 


15,856 


21,56 


6,787 


15,057 


5,227 



(c) Central weight vectors (d) Barycenter 



Students 


Wl 


w 2 


w 3 


Wl 


w 2 


w 3 


S3 


0,61 


0,22 


0,17 


0,25 


0,42 


0,33 


S5 


0,02 


0,42 


0,56 








S7 


0,12 


0,62 


0,26 









Because the Dean has to give the scholarship to only one of these students, and she wants to be as much 
cautious as possible in choosing the best one, she would like to have more information about the possible 
rankings that could be obtained varying the weights assigned to the different subjects. For this reason, we 
decided to show her also the results obtained by applying SMAA to the classical PROMETHEE II method. 
We pointed out that, looking at the final ranking of the students obtained with respect to the net flows, 
student S3 is preferred to student S7 with a frequency of 53.61% while S7 is preferred to S3 with a frequency 



of 43.69%; looking at Table [3(b)] S7 will get always one of the first three positions in the final ranking, while 



S3 will reach the first two positions with a frequency of the 86,65% but she will reach a position between 
the fourth and the seventh with a frequency of 10.13%. 



Sl 
S2 



S 4 
S.5 



S7 

S8 



32,08 




28,74 
6,69 




18,12 
18,68 
31,98 




10,23 
29,64 
23,80 
7,59 




38,92 
51,48 
9,25 
26,25 
41,90 




9,09 
0,80 
25,06 
1,07 
12,39 
15,16 




9,63 
25,24 
27,05 
5,83 


39,59 
10,78 




At the same time, we observed that only students S3, S5 and S7 can get the first position and this happens if 



the most important criterion is Mathematics, Literature or Physics respectively (see Table 3(c) ). Looking at 
Table |3(d)| we get also that independently from which alternative reach the first rank, the most important 
criterion for the DM is Physics followed from Literature and Mathematics. 

In the second part of this didactic example, we shall show how to apply the SMAA methodology to the 
bipolar PROMETHEE methods. 

As the Dean thinks that scientific subjects are more important than humanistic subjects but, at the same 
time, she does not want to favour students good in scientific subjects but having lacks in humanistic sub- 
jects. Besides, she thinks that Mathematics and Physics are redundant because, generally, a student good 
in Mathematics is also good in Physics and viceversa. 

The Dean stated that, locally, student si is preferred to student S2 more than student S3 is preferred 
to student 54, and that student sj is preferred to student ss more than student S5 is preferred to stu- 
dent sq. Translating these information with the constraints Ch B (P B (si, S2), /t) > Ch B (P B (53, S4), fi) and 
Ch B (P B (s-?, sg), fi) > Ch B (P B (S5, sq), £1) and solving the optimization problem (fTTj) . we found that e* < 0, 
therefore the classical PROMETHEE methods are not able to restore the preference information provided 
by the Dean. Solving instead the optimization problem (|12p. we found that e\ > and therefore the bipolar 
PROMETHEE methods are able to restore these preference information. 

Table 4: Results obtained by the bipolar PROMETHEE I method 
(a) Preferences in percentage (b) Incomparability in percentage 





Si 


S2 


S3 


S4 


S5 


se 


S7 


ss 


si 





15,502 








0,636 


6,384 





0,668 


S2 

















2,844 








S3 


64,839 


100 





95,552 


52,853 


99,897 


0,242 


68,591 


S4 


13,584 


86,066 








16,907 


45,775 





30,243 


S5 


74,384 


82,306 





33,646 





26,957 


(1 


51,839 


S6 





11,519 





0,009 














st 


100 


100 


42,288 


100 


99,958 


99,845 





100 


S8 


31,673 


78,226 





30.05 


2,122 


29,272 









Si S 2 S3 S4 S5 S(i S7 S 8 

si 84,498 35,101 86,416 24,98 93,616 67,059 

s 2 13,934 17,694 85,637 21,774 

s 3 4,448 47,147 0,103 57,47 31,409 

s 4 49,447 54,216 39,707 

s 5 73,043 0,042 46,039 

s 6 0,155 70,728 

s 7 

ss 



In Tables 4(a) and 4(b) are shown the results obtained by applying SMAA to the bipolar PROMETHEE 



I method, while Tables 5(a) , |5(b) and Table 5(c) present the results obtained by applying SMAA to the 
bipolar PROMETHEE II method. By using the bipolar PROMETHEE I method, and therefore considering 
separately the bipolar positive flow and the bipolar negative flow respectively, we pointed out again, that 
students S3 and S7 are the best because S3 is preferred to all other students, except to S7, with a frequency of 
at least 52.85% while S7 is preferred to all other students, except to S3, with a frequency of at least 99.84%. 
Comparing S3 and S7, we got that S7 is preferred to S3 with a frequency of 42.28%, S3 is preferred to S7 with 
a frequency of 0.242% and they result incomparable with a frequency of 57.47%. 

By applying the SMAA methodology to the bipolar PROMETHEE II method, we get that considering 



Table 5: Results obtained by the bipolar PROMETHEE II method 



(a) Preferences in percentage (b) First rank acceptability 





si 


S2 


S3 


S4 


S5 


S6 


S7 


ss 




b i 


b ? 


b ? 


b f 


b ? 


b ? 


b i' 


b ? 


si 


n 


100 





10,193 


1,998 


99,267 


n 


2,383 


si 











1,569 


10,938 


87,258 


0,235 





S2 

















32,995 








s 2 




















32,995 


67,005 


S3 


100 


100 





100 


99,803 


100 


11,243 


99,939 


S3 


11,243 


88,56 


0,136 


0,061 














S 4 


89,807 


100 








40,787 


98,613 





50,856 


s 4 








40,161 


11,321 


38,27 


8,916 


1,332 





S5 


98,002 


100 


0,197 


59,213 





100 


n 


86,808 


so 





0,189 


54,347 


36,416 


7,591 


1,457 








se 


0,733 


67,005 





1,387 








n 





S(i 














0,553 


1,014 


65,438 


32,995 


S7 


100 


100 


88,757 


100 


100 


inn 


n 


inn 


S7 


88,757 


11,243 




















S8 


97,617 


100 


0,061 


49,144 


13,192 


100 








S8 





0,008 


5,356 


50,633 


42,648 


1,355 









(c) Central weight vectors 



Students 




a 2 


«3 


ai2 


ai3 


123 


a l|2 


a l|3 


a + 
a 2|l 


s.i 


0,69 


0,62 


0,57 


-0,26 


-0,32 


-0,29 


-0,21 


-0,17 


-0,08 


■17 


0,56 


0,69 


0,64 


-0,29 


-0,26 


-0,34 


-0,14 


-0,16 


-0,11 






a 2|3 


a 3|l 




a l|2 


°1|3 


B 2|l 


a 2|3 


°3|1 


°3|2 




-0,09 


-0,09 


-0,10 


-0,08 


-0,09 


-0,21 


-0,10 


-0,17 


-0,09 


§7 


-0,16 


-0,05 


-0,16 


-0,11 


-0,05 


-0,14 


-0,16 


-0,16 


-0,16 



(d) Barycenter 



a-l 


a-i 


03 


112 


013 


a-23 


a i*|2 


a t\3 


a 2|l 


0,625 


0,655 


0,605 


-0,275 


-0,29 


-0,315 


-0,175 


-0,165 


-0,095 




" 2 |3 


a 3|i 


a t\2 


°1|2 


°1|3 


°2|1 


°2|3 


a 3|i 


°3|2 


-0,125 


-0.07 


-0,13 


-0,095 


-0.07 


-0,175 


-0,13 


-0,165 


-0,125 



the ranking obtained taking into account the bipolar net flow and all sampled sets of parameters, only 
students S3 and S7 can reach the first position. Besides, they shall fill surely the first two positions but s-j 



will be the best with a frequency of 88.75% and, looking at Table 5(c) , just in case Physics is the most 
important criterion, followed by Literature and Mathematics. 

We observe in Table [5(d)] that , for the DM, the criteria are almost equally important and there is redundancy 
between all couples of considered criteria. 



7 Conclusions 

PROMETHEE methods and Stochastic Multiobjective Acceptability Analysis (SMAA) have been widely 
applied to deal with several real world problems [3J, [32] ■ In this paper we have proposed to apply SMAA to 
the classical PROMETHEE methods [4, 5j and to the bipolar PROMETHEE methods [8, 9j. The integration 
of SMAA and PROMETHEE methods permit to get recommendations provided considering several set of 
parameters compatible with the preference information provided by the Decision Maker (DM), and not 
considering only one of these sets. Thus the methodology presented in this paper permits to deal effectively 
with robustness concerns related to the choice of preference parameters in PROMETHEE methods, also in 
case interaction between criteria and bipolarity of the scales are considered. 

In the case of PROMETHEE II and of the bipolar PROMETHEE II, one can look at the possible final 
rankings obtained considering the net flow and the bipolar net flow respectively and at which parameters 
give to an alternative the best position. 

In the didactic example we have shown that the SMAA methodology can be applied to PROMETHEE 



methods, also in cases the criteria are interacting. This brings us to believe that the methodology we are 
proposing in this paper can be adopted to deal with many complex real-world problems. 
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